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Abstract. A nonequilibrium liquid theory for uniformly sheared granular liquids is developed starting from SLLOD equa- 
tions. We derive a generalized Green-Kubo formula and demonstrate that it yields the nonequilibrium steady-state average 
which is essentially independent of the choice of the initial condition. It is also shown that the fluctuating hydrodynam- 
ics can be derived from Mori-type equations for density and current-density fluctuations if one considers a weak-shear and 
small-dissipation limit along with the Markovian approximation. 
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INTRODUCTION 

Assemblies of dense granular materials behave unlike usual materials. One may think that a liquid phase is absent 
for granular assemblies since they lack attractive interaction. Indeed, the gas kinetic theory such as the Boltzmann- 
Enskog theory has historically been employed in describing dense granular flows, where correlation effects appear 
only through the contact value of the radial distribution functional, 2, 3, 4, 5, 6] This approach is powerful and semi- 
quantitatively accurate even for considerably dense systems. In these days, however, we have recognized the relevancy 
of the concept of "granular liquids" because correlation effects in granular flows, such as long-time correlations and 
long-range correlations, turned out to be relevant as in molecular liquids. [7, 8, 9, 10, 1 1, 12, 13] 

There have been some developments in the granular liquid theory starting from microscopic basic equations such as 
the Liouville equation. The Liouville equation of granular fluids was first discussed by Schofield and Oppenheim[14] 
long time ago, and Brey et al.[\5] further developed such a formulation. Recently, Dufty et al.[16, 17] discussed in 
detail the Green-Kubo formula of freely cooling granular gases starting from the Liouville equation. The Green-Kubo 
formula for granular fluids has been discussed in various contexts[18, 19, 20, 21, 22], and these studies suggest that 
some correction terms to the conventional Green-Kubo formula are necessary for granular fluids. 

There is an advantage in using the liquid theory. It is known that liquid theories such as the mode-coupling theory 
(MCT) for supercooled liquids are commonly used to describe the glass transition of molecular liquids or colloidal 
assemblies. [23, 24] Liu and Nagel[25] proposed that the jamming transition is a fundamental transition in glassy and 
granular materials. Since then many aspects of similarities between the conventional glass transition and the jamming 
transition have been investigated [26], where some researchers have used granular materials to study dynamical 
heterogeneity in glassy materials. [27, 28, 29, 30, 31] Along this line it appears natural to apply a liquid theory to 
dense granular assemblies. 

Nevertheless, the jamming transition under a plane shear exhibits some distinct aspects from conventional glass 
transitions, e.g., the jamming transition depends strongly on details of microscopic interactions between particles. In 
fact, one cannot use a naive MCT [32] in describing sheared jamming transitions. Instead, the jamming transition of 
frictionless granular particles is believed to be a continuous transition at a critical density above which elastic moduli 
and the yield stress become nonzero, and there are scaling laws in the vicinity of the critical point as observed in 
conventional critical phenomena. [33, 34, 35, 36, 37, 38] 

Quite recently, Chong and Kim[39] have reformulated the liquid theory of sheared dense molecular liquids with the 
Gaussian thermostat. Later, Chong et al.[40] extended their formulation to soft granular liquids under a plane shear, 
and found the existence of the generalized Green-Kubo formula and integral fluctuation theorem. Chong et al.[41] have 
further developed an MCT for uniformly sheared granular liquids. This generalized Green-Kubo formula applies not 
only to linear regime but also to nonequilibrium states arbitrarily far from equilibrium. Thus, one can avoid the use of 
hydrodynamic equations including Burnett and super-Burnett terms which occasionally exhibit a divergent behavior. 



The aim of this paper is to outline our formulation of the liquid theory for dense sheared granular materials. In 
the next section, we summarize the microscopic starting equations, such as SLLOD equations, Liouville equations, 
and nonequilibrium distribution function. In section III, we present the transient time correlation function formalism 
and discuss steady state properties. We then derive the generalized Green-Kubo formula. It is demonstrated that the 
generalized Green-Kubo formula yields the nonequilibrium steady-state average which is essentially independent of 
the specific choice of the initial condition. In section IV, we outline the liquid theory beyond the generalized Green- 
Kubo formula, and discuss its connection to the fluctuating hydrodynamics. The paper is summarized in section V. 



MICROSCOPIC STARTING EQUATIONS 

In this section, we derive exact microscopic equations and relations which serve a basis in constructing a nonequilib- 
rium liquid theory for uniformly sheared frictionless granular particles. 



SLLOD equations of motion 

Let us consider a system of N smooth granular particles of mass m in a volume V under a stationary shear 
characterized by the shear-rate tensor k. We assume that each granular particle is a soft-sphere, and the contact 
force acts only on the normal direction. Under a homogeneous shear, the velocity profile is given by k ■ r at position 
r. The Newtonian equations of motion describing such a homogeneously sheared system are given by the SLLOD 
equations [42] 

t i = ^- + K-r i , (1) 
m 

where r, refers to the position of the ;'th particle, r, = dVi/dt, and 

p i =Ff»+F|™'- 1 f.p i . (2) 

(el) (el) 

Here Eq.(l) is the definition of the peculiar momentum p, satisfying £ ( p, = 0, F- = L^ F^. is the conservative 
force exerted on the ;'th particle by other particles with 

F (eD = _ a« = e(ff rik)f[d ^ (3) 

where a is the diameter of each grain, t<(r (J t) is the pairwise potential, r,* = r, — r^, = |r ( -^|, = r^/r,^, and &(x) 
is the Heviside function satisfying ©(x) = 1 forx > and &(x) = forx < 0. The actual elastic repulsive force f{x) 
is proportional to x 3//2 for three dimensional systems, but we sometimes use a simpler form f(x) °c x. Similarly, the 
viscous dissipative force f| v1s ' 1 is represented by a sum of two-body contact forces as f| v1s ' = L 7 y,Fjj ls ' with 

F|f ' = -*iJ*(rij)(BiJ-*ij) = -iiMd-rtjMo-rijHgij-ttj). (4) 

In Eq. (4) we have introduced gy = v,- — v,- = (p,- — Pj)/m + K ■ (r\ — rj ) with the velocity of ith particle v; = c/r,-/ dt . 
The viscous function y(x) is proportional to y/x for three dimensional systems, but we sometimes use a simpler model 
y(x) = const. 

It should be noted that SLLOD equations (1) and (2) reduce to the Newtonian equation of motion 

mr i =Ff+F\ vis \ (5) 

if one eliminates the peculiar momentum p,. This means that the SLLOD equations are equivalent to the Newtonian 
equation of motion under the Lees-Edwards boundary condition. [42] We also note that frictionless sheared granular 
particles under a constant pressure boundary behave as those under the Lees-Edwards boundary condition in the 
vicinity the jamming transition. [43] On the other hand, it is hard to extract physical essences from actual frictional 
granular assemblies under a physical boundary condition. Thus, a set of SLLOD equations is a natural starting point 
in constructing a liquid theory of granular particles. We also notice that a granular flow on an inclined slope can be 
approximately described by a uniformly sheared flow except for the boundary layers. [44] 



We stress that an energy sink term is necessary even for a system of a sheared molecular liquid. Indeed, the system 
heats up without the energy sink. In simulations, one usually introduces a thermostat, and an experimental apparatus 
for the real system plays a role of the thermostat. We have already confirmed that behaviors of uniformly sheared 
granular liquids with small inelasticity are almost the same as those for a model of a molecular liquid with a velocity 
rescaling thermostat. [1 1, 12] 

We have assumed that the interaction between granular particles is described by a soft-core model, which differs 
from most of the conventional treatment of the granular gas kinetic theory and the MCT for sheared granular liquids 
where a hard-core model is adopted. [32] The soft-core model has several advantages: (i) the model is more realistic 
than the hard-core model, (ii) one can apply this formulation to very dense systems even in the vicinity of the jamming 
transition, and (iii) one can avoid the use of the pseudo-Liouvillian. 



The Liouville equation 

For nonequilibrium systems described by the SLLOD equations, the Liouville equation is commonly used. [42] The 
time evolution of phase variables whose time dependence comes solely from that of the phase space point T = (r^, p N ) 
is determined by 

^A(T)=f~A(r) = iJ2fA(r). (6) 
The operator i££ is referred to as the Liouvillian. The formal solution to this equation can be written as 

A(rv) = exp(m)A(r). (7) 

On the other hand, the Liouville equation for the nonequilibrium phase-space distribution function p (T, t) is given 

by 



dp(T,t) 



dt 

where the phase space contraction factor A(T) is defined by 

d ^ «-,/ d 



t~+A(T)]p(T,t) = -i^p(T,t), (8) 



(9) 



For our model (1) - (4) one easily obtains its explicit form: 



A(r) = — -£©(«/- ry)y(£/-r y ) < 0. (10) 

m U 

The formal solution to the Liouville equation (8) reads 

p(r,f)-ex P (-/^ t f)p(r,o). (ii) 

From Eqs. (6) and (8) we readily obtain the relation 

iJSf t (T) = iJSf(r)+A(r). (12) 
One can show that the following adjoint relations hold [42]: 

J dY[i^A{Y)]B{Y) = -J dFA(r) [;Jz? f B(r)], (13) 

yrfr[e m A(r))]B(r) = JdrA{r) e - iJft 'B(r). (14) 

If the phase-space contraction factor A(r) is identically zero, then i££ ' = ;'Jz? holds, and the Liouvillian becomes 
self-adjoint or Hermitian. In general, this is not the case for nonequilibrium dissipative systems. 



Nonequilibrium distribution function 



Let us consider an equilibrium system to which a constant shear rate y satisfying K a p = ySaxSpy is applied at time 
t = 0, and thereafter the system evolves according to the SLLOD equations (1) and (2). The Liouvillian is given by 
iJgf = ;J^ (el) for t < and 

iSe = + iSff + iJ? {v[s) (15) 

for t > 0, where the unperturbed adiabatic or the elastic part (i_Sf( e1 )), the shear part (i-Sfy), and the viscous part 
(i'Jz? ( V1S )) are respectively given by 



Hi. A +F (el) .— 

m <?r; 1 dp,- 



iJSf* 



■^(vis) 



(vis) 



dp,' 



Here, we assume that the initial distribution is given by the canonical one 



p(r,o) = Peq (r) 



Z(J3) = JdTe-W, 



Z(P) ' 

where /3 = 1/T is the inverse temperature in the initial state, and H is the total Hamiltonian defined by 



"=L< 



2m 2 



(16) 
(17) 
(18) 

(19) 

(20) 



It should be noted that the effect of the shear appears in H through pf /2m = m(v, — re- r;) 2 /2. We also note that there 
holds a trivial relationship 

;^ (el) Pe q = 0. (21) 

It might appear that our formulation depends strongly on our choice of the initial condition, Eq. (19). However, we 
will argue below that nonequilibrium steady-state properties therefrom are insensitive to such a choice. 
From Eq. (8) the time evolution of the distribution function can be written as 



With the identity 

Eq. (22) can be expressed as 



-LSfh 



p(r,0=e- Urt, peq(T). 

1+ f dse-^'l-itf), 
Jo 



p(r,o = Peq(r)+^^ e -^(-^t )Peq(r) . 

From Eqs. (10), (12), (15)-(18), and (21) we get 

^ t Peq(r) = /^pe q (r) + /^( vis ) Peq (r) +A(r) Peq (r). 

The first term on the right hand side in this expression is given by 

^ Peq (r) = /3£[(K.r i ).F i (el) + (K. Pi ).Hi]pe q (r) 

= rc:cr( el )pe q (r) = yai y el )pe q (r) 



(22) 



(23) 



(24) 



(25) 



(26) 



where a' e1 ' denotes the elastic stress tensor whose element is given by 



< ) =D^+^ (el)/3 ] («>P =*,**)■ (27) 



In the final equality of Eq. (26) we have used the specific form K a p = y8 ax 8^ y . On the other hand, from Eq. (18), 
/^ vis >p eq (r) is given by 

/j^Wr) = -/3p eq (r)£F,^.Pi 

i m 

= p to^-Pi) -F^Wr) - § £F,f ) .g ikPeq (r). (28) 

i Z i,k 

It is convenient to introduce Rayleigh's dissipation function Si as 



(29) 



p(r,/) can then be written as 

P(T,0 - Peq(r) + ^^-^ ts [p eq (r)a(r)], (30) 

where £2(T) is the nonequilibrium work function defined by 

n(r) = -/3rcj^(r)-2/3^(r)-A(r) (3i) 

in terms of the total stress tensor 

° a p=l} P -^+r?{F^+F^)] (a,P=x,y,z). (32) 

i m 

One can show that the equilibrium average of the nonequilibrium work function is zero: 

(Q(T)) c =0. (33) 

Hereafter, we shall reserve the notation (• • - ) c for representing the averaging over the initial canonical distribution 
function p(T,0) =peq(r): 

(■■■)c = JdTp cq (T)---. (34) 



GENERALIZED GREEN-KUBO FORMULA 
Transient time-correlation function formalism and generalized Green-Kubo formula 

In contrast to the case of equilibrium quantities, the nonequilibrium ensemble average (A(t)) c of a phase variable 
A(t) depends explicitly on the time t past since the start of the shearing. Using the nonequilibrium phase-space 
distribution function p(r,f), (A(r)) c can be expressed as 

(A(f)) c = j dTp{T,Q)A{t) = J dTp{T,t)A{0), (35) 

where the second equality follows from Eq. (14). 

Substituting Eq. (30) into Eq. (35) and then using Eq. (13), one obtains 

(A(0)c = (A(0))c+ f ds(A(sM0)) c . (36) 
Jo 



The expression (36) relates the nonequilibrium value of a phase variable A (t) at time t to the integral of the transient 
time-correlation function. Indeed the integrand in Eq. (36), (A(s)Q.(0)) c , is the correlation between the nonequilibrium 
work function in the initial state and A at time s after the shearing force is turned on. It should be remembered, however, 
that the dynamics inside the brackets (• • -) c is governed by the granular SLLOD equations, and only averages like 
(A(0)) c coincide with equilibrium quantities. 

The system is in a nonequilibrium steady state if the ensemble averages of all phase variables become time- 
independent. Let us notice that the long-time limit of Eq. (36) approaches a constant, and hence, the integral is 
convergent if the system displays mixing. [42] This feature can be demonstrated by taking a time derivative of Eq. (36): 

^<A(0)c = <A(0n(0))c (37) 

If the system displays mixing [42], all long-time correlations between phase variables vanish. With the aid of Eq.(33) 
we obtain \m\ t ^ oa {d/dt)(A{t)) c = and 

lim(A(f)) c = (A) ss . (38) 

t— 5-°° 

Here, the steady-state average, denoted by (• • -) ss hereafter, is obtained from the long-time limit of Eq. (36): 

(A) ss = (A(0)) c + H ds(A(s)£l(0)) c . (39) 
Jo 

Because of Eq.(33), Eq.(39) can be rewritten as 

(A) ss = (A(0)) c + f ds(AA{s)a{0)) e , (40) 
Jo 

where AA(s) = A(s) — A(s -> °°). 

Equation (39) is the generalized Green-Kubo formula which relates the steady-state average to the time-correlation 
function describing transient dynamics evolving from an initial equilibrium towards a final steady state. One can easily 
show that Eq. (39) reduces to the conventional Green-Kubo formula if the external force is weak and the dissipative 
force is neglected, i.e., for small y and y(x) — 0. For example, by setting A = oVy in Eq. (39), one obtains for the 
steady-state shear stress defined via a ss = —(<Jxy)ss/V 



K(o))c 



-±jH<fa<a^(j)Q(0)) c . (41) 



V 

When y(x) = 0, there hold (<7xy(0)) c = and Q. = -fijOxy [see Eq. (31)], and Eq. (41) formally reduces to 



%f f dj<a^(j)a^(0)) c . (42) 
V Jo 



For small y, one can replace the Liouvillian governing the dynamics of a v (s) in the integrand by that for a quiescent 
equilibrium state, and hence, Eq. (42) is the conventional Green-Kubo formula for the viscosity r) defined via 
tj = <7 ss /y. It should be noted that the conventional derivation of Green-Kubo formula requires a convergent factor 
e -et f or tne integrand with taking the limit e — > 0. Similarly, a small dissipation is also necessary for Eq.(41) to obtain 
a convergent result. Namely, the dissipation plays a role of the convergent factor. Otherwise, the system is heated up 
and cannot reach a steady state. 

Therefore, with the aid of Eq.(41), the viscosity 7] satisfies 

t l = -&A-^£* {<hymme (43) 

in general situations. We should stress that Eq. (41) or Eq.(43) is the full-order expression, and applies to nonequilib- 
rium states arbitrarily far from equilibrium. Thus, we do not have to worry about Burnett or super-Burnett terms which 
occasionally exhibit an unstable behavior. In other words, the viscosity 77 or the steady shear stress (7 SS involves effects 
of nonlinear rheology, and is free from the magnitude of the deviation from a reference state. 



On the initial condition 



In this subsection, let us demonstrate that (A) ss is independent of the choice of the initial condition such as the 
initial temperature and initial distribution. [40] This result is highly nontrivial, because Eq. (39) appears to depend on 
the choice of the initial canonical distribution. 

From Eqs. (1 1) and (19) one obtains the Kawasaki representation [42] 



(44) 



r [' 1 e-^H) r ft 

p(r,f) ^expj-^ dsA(-s)\ =peq(r)exp^ ds£l(-s) 

In the second equality we have introduced the nonequilibrium work function Cl(t) = e'-^'Cl(r) at time t. 
Since 

^{^w^ m - H{ - t)] ^w (45) 

one obtains from Eqs. (44) and the definition of the average for t — > °° 

^(A(0)c = JdTA(f))[{H) e -H(-t)]p<r,t) 

= (A(t)) c (H) c -(A(t)H) c ^0, (46) 

i.e., (A) ss = lim r ^oo(A(f )) c is independent of the inverse temperature j3 of the initial equilibrium state. Thus, (A) ss is 
uniquely specified by the "thermodynamic" parameters (N,V, y) characterizing the nonequilibrium steady state. 

One can prove a stronger statement that the average of any variable in the steady state is invariant if the initial 
condition can be expanded in an orthogonal polynomial of the kinetic energy associated with the Gaussian function or 
the exponential function such as the Laguerre bi-polynomial and the Hermite polynomial. Let us demonstrate that the 
average under the initial condition expanded in the Sonine polynomial, which is related to the Laguerre bi-polynomial, 
is the same as the one under the canonical initial condition. We assume the following initial condition 

p..(n^{i + |< 2 (o£)}. (47) 

where «/ and Sy 2 (x) are respectively the expansion parameter and the Sonine polynomial which satisfies the orthogo- 
nal condition 



L 



^ dxe- x x m S%>(x)S%>(x)= v ' 8 M (48) 

with the Gamma function T(x). It should be noted that the Sonine expansion around the Gaussian has widely been used 
for the description of freely cooling granular gases[2], but Eq.(47) is more general than the case of freely cooling cases. 
Indeed, any function of the kinetic energy can be expanded in an orthogonal polynomial. Thus, the only assumption 
adopted here is that the initial condition can be represented by a product of the canonical distribution and a function 
of the kinetic energy. 

Let us denote the average under Pi n (r) as 

(A(f)) in = J dTp ia (T)A(t). (49) 

We also introduce the difference between the average under the initial condition (47) and the one under the canonical 
initial condition (19): 

5A(f) = <A(f)>in-<A(f)>c (50) 



One immediately obtains 



SA(t) = |jrA(r)( Pin (r(0)-pc(r(0)) 

= |jrA(r,0(Pin(r)-Pc(r)), (5i) 
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FIGURE 1. The time evolution of the granular temperature starting from To = 5 x 10 6 , Tq = 5 x 10 5 , and the homogeneous 
cooling state. See the text in details. 



where we have used Eq.(14). Substituting (47) into (51) we obtain 

= |„,^,r,,) S ™ ( ? |l))^| ,^,r,,» t (5«(^)) t (52) 

in the limit t — > °° with the aid of the mixing property. If we use Sy 2 (x) = 1 and Eq. (48), we get the relation 




Thus, we obtain 

lim<5A(f)=0, (54) 

t— p- 00 

which is the end of proof. Thus, the steady state starting from Eq.(47) is equivalent to the one from the canonical 
distribution (19). 

To demonstrate the irrelevancy of the choice of a specific initial condition, we shall present computer-simulation 
results for two-dimensional soft granular particles. The simulations have been done for the canonical initial condition 
and for the homogeneous cooling state. The system consists of polydisperse 4000 grains of diameters 0.1 Oq, 0.8on, 
0.9ob and on> and the number of grains of each diameter is 1000. The total area fraction is = 0.5. We have adopted 
the linear spring model to represent the elastic repulsion during a contact. All variables are non-dimensionalized by 
the maximum diameter of grains Co, its mass m and the spring constant k. The dissipation appears through the linear 
viscous damping with its coefficient Tj^amp = 1-0. (m, k, and 77damp = 10 are common to all the grains.) The applied 
shear rate is y = 0.0005. The freely cooling initial condition has been prepared by performing a granular simulation in 
the absence of shear up to t = 200 starting from the canonical distribution at Tq = 5 x 10~ 5 . 

Figure 1 shows the result of the granular temperature defined by the kinetic energy. It is easily seen that all the 
results starting from different initial conditions converge to a unique steady kinetic temperature in the long time limit. 

Figure 2 displays the radial distribution function g(r) for the largest grains in a steady state. We find that the radial 
distribution function g{r) in the steady state is independent of the initial condition. Thus, our numerical results support 
our theoretical prediction that steady-state properties are independent of the choice of a specific initial condition. 

Thus, the generalized Green-Kubo formula (41) yields the same steady-state average irrespective of the initial 
condition. We should note that most of nonequilibrium generalizations of the Green-Kubo formula assume the 
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FIGURE 2. The radial distribution function in the steady state. The legend is common with that in Fig.l. 



existence of a nonequilibrium steady distribution function and strongly depends on its steady distribution p ss . This 
standard method has several difficulties such as (i) the determination of p ss is difficult, and (ii) the distribution 
might not approach a steady value, because there are in general no compatible solutions of both ;JzfA ss (T) = and 
/Jz? T p ss (r) = 0. On the other hand, our formulation is free from such difficulties, and we can calculate, e.g., the steady 
shear stress under the canonical initial condition and the obtained result is independent of the adopted initial condition. 



OUTLINE OF THE LIQUID THEORY BEYOND GREEN-KUBO FORMULA 

In this section, we briefly explain how to use the generalized Green-Kubo formula to describe sheared granular liquids. 
Because of the limitation of the length of this paper, we shall skip details of the derivation, and the interested reader is 
referred to ref.[41]. 

Because the energy is not a conserved quantity which quickly relaxes to a steady value in the uniform shear, the 
relevant hydrodynamic variables are the density fluctuations 

n q (0 = £A r < W -ATS q , , (55) 
and the current density fluctuations defined by 

4=I^ q - ri - (56) 
We introduce the projection operator & onto these variables: 

where Sk = («k( f )«k( f )*) /N and vj = y/T/m. The complementary projection operator shall be defined by =3 = /— S?. 
Let us introduce 

R\{t) = e^ at R\ (58) 

with 

2 

R\ = m£j\=i£j\-iq x V ^n < ,-iB\n ix +A\»j^ (59) 



where i& = ;Jz^ el ) + L5f < vis >, and 
m J 

iB\ = --^-{n J drf x xyr,^{r)g{r)e iCiI +n 2 J dr J dr' r x xyr^(r)g^ (r^e 1 ^^, 
in terms of the pair and triple correlation functions 



ng(r) 



^L<*('-'«)>. 



«V 3) (r,r') = - £ (5(r-r y )5(r'-r a )>. 



(60) 
(61) 

(62) 
(63) 



One finds the following continuity equations for the sheared system relating the partial time derivative of n q (t) = 



• iJft n q to 4(0=«""^: 



dr ^ <?q 



«q(0 = iq-j q (0, 



(64) 



and 



5 <? 1 

— -q- K- — 
<9r M <?q 



4(0 



^n q (r) + a^«q(0 -A^m+R\{t) fidsMptfeW-'^ 
+ f dsiL\ {s) e'^'-\ {s) f dsN^ (s) e>*l'-% s) , 



where we have introduced the following memory kernels 



K{t) = 



-''jv^7^(O^K w "(o)]>. 



(65) 

(66) 
(67) 
(68) 



The equations (64)-(68) are the exact equations for uniformly sheared granular liquids governed by Eqs. (l)-(4). 

To obtain a closure to these equations, one has to introduce some approximations such as the mode-coupling 
approximation. The details of such approximations will be reported elsewhere. [41] Instead, here, let us briefly explain 
what equations can be obtained under the Markovian approximations in the case of a weak shear and an elastic limit. 

It is straightforward to show that Eq.(65) reduces to 



— q k — 

dt dq 



4(0 « iqxfn q ( t )-£dsM^(s)e'^f q(s) +R X (t) 



(69) 



under the weak shear and the elastic limit. Thus, the effect of shear appears only through the convective deformation 
of the wave number in the elastic and the unsheared limit. This equation still includes the non-Markovian memory 
kernel. 

For many situations in a liquid state far from the jamming transition, one can ignore memory effects. When we adopt 
the Markovian approximation with the assumption that the system is isotropic, it is known that the memory kernel can 
be approximately given by its hydrodynamic limit[45] 



(70) 



where qx = qx /q, V\ and V2 are the bulk kinetic viscosity and the kinetic viscosity, respectively. Equation (69) then 
reduces to 



c) d - . v? 

S, 



dt dq 



jq(t)^iqx^n q (t)-V 1 qx(q-j (l (t))~V2q 2 (8x^-qxqn)jq(t)+Ri l (t), (71) 

q 



where R q (f ) satisfies the fluctuation-dissipation relation 

(<(f)<(0)) ^8(t)8(q + k)Nv 2 T {qxq^q 2 Vl + (8x^-qxq^q 2 V2}. (72) 

Equation (71) combined with Eqs. (64) and (72) is the equation of fluctuating hydrodynamics. 

Let us notice that both the equal-time long-range correlation function[12] and the long-time tails of autocorrelation 
functions[ 1 1,13] can be discussed within the framework of the fluctuating hydrodynamics. Our liquid theory presented 
here, therefore, provides not only a microscopic basis of the fluctuating hydrodynamics, but also a basis of both the 
long-time tails and the long-range correlations for sheared granular liquids. 



DISCUSSION AND CONCLUSION 
Discussion 

Now, let us compare our formulation with previous formulations of the Green-Kubo formula for granular fluids. [16, 
18, 19, 20, 21, 22] We first stress that our formulation is unique in that it applies to a nonequilibrium steady state of 
uniformly sheared granular liquids, while the previous ones deal with the Green-Kubo formula for granular gases in 
a freely cooling state. We also note that Green-Kubo formula in the previous studies is the linear response theory to a 
nonequilibrium "steady state" (which includes a homogeneous cooling state), but our method is a nonlinear response 
theory to the initial canonical state. In practice, the application of the linear response theory to freely cooling granular 
gases has some technical problems; (i) one cannot take the t — > °° limit of the integral of the time-correlation function 
because the granular particles quickly lose their kinetic energy, (ii) thus the behavior of the steady state strongly 
depends on the cut-off time of the integration, and (iii) the determination of a nonequilibrium steady distribution is 
difficult. Let us also notice that that freely cooling states cannot be realized in experiments. On the other hand, our 
method which is the nonlinear response theory to the initial canonical distribution has several advantages such as (i) 
the steady state under a uniform shear can be approximately realized in many situations, and (ii) the steady state is 
almost independent of the choice of the initial condition, although our method cannot be generalized to a response 
theory to a reference state. 

Nevertheless, it is remarkable that the Green-Kubo formula for freely cooling granular gases has a similar structure 
to our generalized Green-Kubo formula. Indeed, Eq. (139) in ref.[16] is the essentially same as Eq.(39) where their 
expression corresponds to the integral Tj °c J^dt (Oxy(0)[e~ l ^'Q.(r)]} in our context. 



Conclusion 

This paper summarizes our recent studies on the granular liquid theory under the uniform shear. We demonstrated 
that there exists the generalized Green-Kubo formula in sheared granular liquids. We also showed that it yields the 
nonequilibrium steady-state average which is essentially independent of a specific choice of the initial condition. It 
should be noted that our formulation does not rely on the presence of the steady distribution function which is hard 
to obtain, but any averaged quantity in the steady state can be obtained under the simple initial condition. In the 
previous section, we outlined how we can apply our formulation to characterize the behaviors of granular liquids. One 
of the most important conclusions in this paper is that the granular liquid theory starting from the Liouville equation 
reduces to the fluctuating hydrodynamics if one considers a weak-shear and elastic limit along with the Markovian 
approximation. Thus, this paper provides a microscopic support of the fluctuating hydrodynamics which is known to 
give accurate results of the time correlations and the spatial correlations. 

We have skipped detailed derivation of Mori-type generalized Langevin equation and MCT This will be discussed 
in another paper. [41] We also note that we can derive the integral fluctuation theorem along the same line, though there 
is no microscopic time reversal symmetry. [40] The formulation of granular liquid theory, thus, gives an interesting 
subject in pure nonequilibrium statistical mechanics. [46] 
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